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Abstract 

We have carried out a high statistics (2 x 10 9 events) search for ultra-high energy gamma-ray 
emission from the X-ray binary sources Cygnus X-3 and Hercules X-l. Using data taken with the 
CASA-MIA detector over a five year period (1990-1995), we find no evidence for steady emission 
from either source. The derived 90% c.l. upper limit to the steady integral flux of gamma-rays 
from Cygnus X-3 is $(£ > 115 TeV) < 6.3 x 10 -15 photons cm 2 sec , and from Hercules X-l it 
is $(£' > 115 TeV) < 8.5 x 10 -15 photons cm -2 sec -1 . These limits are more than two orders of 
magnitude lower than earlier claimed detections and are better than recent experiments operating 
in the same energy range. We have also searched for transient emission on time periods of one day 
and 0.5 hr and find no evidence for such emission from either source. The typical daily limit on the 
integral gamma-ray flux from Cygnus X-3 or Hercules X-l is ^daiiy^ > 115 TeV) < 2.0 x 10 -13 
photons cm -2 sec -1 . For Cygnus X-3, we see no evidence for emission correlated with the 4.8 hr 
X-ray periodicity or with the occurrence of large radio flares. Unless one postulates that these 
sources were very active earlier and are now dormant, the limits presented here put into question 
the earlier results, and highlight the difficulties that possible future experiments will have in 
detecting gamma-ray signals at ultra-high energies. 

1 Introduction 

Cosmic ray particles span a remarkable range of energies, from the MeV scale to more than 10 20 eV 
(eV = electron Volt). At energies above 1 TeV (10 12 eV), we know that cosmic rays do not originate 
from local sources in or nearby our Solar System. Therefore, high energy cosmic rays must come 
from acceleration sites in the Galaxy at large or from outside the Galaxy. Remarkably, after many 
years of research, the exact sites of high energy cosmic ray acceleration remain unknown. 

There are several difficulties that plague efforts to pinpoint the origins of high energy cosmic 
rays. First, since the bulk of the cosmic radiation is electrically charged, any source information 
contained in the directions of the arriving particles is lost due to deflection in the Galactic mag- 
netic field. A second difficulty concerns the energetics of the proposed cosmic ray acceleration 



mechanisms. For example, although models based on shock acceleration in supernova remnants 
offer plausible explanations for the cosmic ray origin up to 10 eV (and perhaps up to 10 15 eV), 
these models become less satisfying and less realistic at energies above 10 15 eV. Since cosmic ray 
origins remain mysterious, it is natural to search for neutral radiation from point sources which, if 
seen, could pinpoint possible acceleration sites. The question of cosmic ray origin is thus a prime 
motivation for high energy neutral particle (gamma-ray or neutrino) astronomy. 

In addition to supernova remnants, possible galactic sources of high energy particles include 
pulsars and compact binary systems. Gamma-radiation has been unambiguously detected from 
the Crab Nebula (a supernova remnant) at energies up to lOTeV by ground-based detectors 

but, historically, the compact binary sources Cygnus X-3 and Hercules X-l have received 
considerably greater attention in the ground-based astronomical community. In the period 1975- 
1990, literally dozens of gamma-ray detections from Cygnus X-3 and Hercules X-l were reported 
by numerous experiments. The detections spanned a wide range of energies (100 MeV to 10 17 eV) 
P], were generally of low statistical significance (typically three to four standard deviations), and 
episodic in nature. Often a statistically significant signal could only be extracted as a result of a 
periodicity analysis, where the data were phase-locked to a known source X-ray periodicity. In spite 
of these difficulties, the sheer number of reports made it difficult to dismiss the detections as being 
entirely due to statistical fluctuation j|. In fact, by the late 1980's, it was generally established 
that Cygnus X-3 and Hercules X-l were powerful emitters of high energy gamma-rays (although 
contrary interpretations of the data existed ) . A number of new, more sensitive ground-based air 
shower arrays were commissioned at this time, including the CASA-MIA experiment in Dugway, 
Utah (USA). This paper describes long-term (1990-1995) observations of Cygnus X-3 and Hercules 
X-l by CASA-MIA. 

In the following section, we outline the experimental techniques of gamma-ray astronomy. 
We summarize the properties of Cygnus X-3 and Hercules X-l, and review previous observations 
of these sources and the astrophysical ramifications from the observations. We then turn our 
attention to the experimental apparatus, the event reconstruction procedures, and the methods 
used to select gamma-ray candidate events. We present results from a data sample of 2 x 10 9 
events, and conclude with comparisons of our results to those of earlier and contemporaneous 
experiments. 

2 Experimental Techniques of Gamma-Ray Astronomy 

Gamma-ray sources typically exhibit power law spectra; fluxes fall rapidly with increasing energy. 
Space-borne experiments (on satellites or balloons) currently have sufficient sensitivity to detect 
gamma-rays up to an energy of ~ 10 GeV. Astronomy at higher energies requires very large 
collection areas available only to ground-based telescopes. Ground-based instruments rely upon 
the fact that high energy gamma-rays interact in the Earth's atmosphere to produce extensive 
air showers. At energies near 1 TeV, atmospheric Cherenkov telescopes use optical techniques 
to detect the Cherenkov radiation in the shower. At higher energies (~ 10 TeV and above), the 
charged particles in the shower penetrate to ground level. Here, air shower arrays record the 
arrival times and particle densities of the charged particles. At energies above 10 17 eV, there is 
enough energy in the shower to allow the detection of nitrogen fluorescence in the atmosphere. 
The faint near ultraviolet fluorescence signal can be optically detected at night by experiments 
such as the Fly's Eye. 
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3 Discussion of the Sources and Earlier Results 



3.1 Cygnus X-3 

Cygnus X-3 is one of the most luminous X-ray sources in our Galaxy ||. The X-ray emission is 
characterized by a 4.8 hr periodicity, which is assumed to be associated with the orbital motion 
of a compact object (neutron star or black hole) around its binary companion. The periodicity 
has been well studied; a complete ephemeris is available for the period from 1970-1995 pL |7], ||, |{|. 
In addition to being a powerful X-ray source, Cygnus X-3 is seen in the infrared and is a strong 
and variable radio source. Radio flares have been detected in which the output from the source 



increases by two to three orders of magnitude on the time scale of days |10|, |ll|, |12|, |l3j. These 
flares were first detected in 1972 and the outbursts have continued through 1994. Since Cygnus 
X-3 lies in the galactic plane, its optical emission is largely obscured by interstellar material. The 
lack of a strong optical signal makes determination of the distance to Cygnus X-3 difficult, but 
general consensus places it near lOkpc [14]. 



The first published result claiming the detection of gamma-ray emission from Cygnus X-3 
came in 1977 from the SAS-2 satellite at low gamma-ray energies (E> 35MeV) |0|. This result 
made use of an apparent correlation between the gamma-ray arrival times and the 4.8 hr X-ray 



periodicity. Later observations by the COS-B satellite [16| with a much larger source exposure 
failed to confirm the SAS-2 result, and the COS-B authors argued that the initial detection [15] was 
flawed because of an incorrect treatment of the diffuse gamma-ray component |jl7fl . To complicate 
matters, there have been re-examinations of both the SAS-2 [[j^] and COS-B |19([ data sets which 
claim that signals indeed exist in both cases. Most recently, the EGRET experiment on the 
Compton Gamma-Ray Observatory (CGRO) failed to detect gamma-ray emission from Cygnus 
X-3 at a sensitivity level comparable to COS-B p0| . To summarize, there exists some controversy 
as to whether low energy gamma-rays have ever been detected from Cygnus X-3. Regardless, it 
can be reasonably concluded that the source is not a strong emitter of gamma-rays in the energy 
range between 30 MeV and 10 GeV. 

The first published report of very-high energy gamma-ray emission from Cygnus X-3 came 
from the Crimean Observatory using an atmospheric Cherenkov telescope at energies above 2 TeV 



[21]. This result was based on data taken between 1972 and 1977 and the emission was claimed 



to be correlated with the 4.8 hr X-ray period. From 1980 to 1990, there were numerous additional 



detections of Cygnus X-3 by atmospheric Cherenkov telescopes [|22j, |23|, |24j, |25|, [26], |27J, |28j. The 
detections were generally episodic in nature and usually required the use of the 4.8 hr periodicity 
to extract a signal. Evidence for a 12.6 msec gamma-ray pulsar inside the Cygnus X-3 system was 



claimed on more than one occasion |2q, 28 



At the higher energies accessible by ground arrays, evidence for ultra-high energy gamma-ray 
emission from Cygnus X-3 was presented by the Kiel array |3(| and subsequently by the Haverah 
Park experiment |il|. These results were based on data taken between 1976 and 1980. The 
gamma-ray emission was apparently steady over this time period and was correlated with the X- 
ray periodicity. Additional evidence for gamma-ray emission from Cygnus X-3 was later reported 



by other air shower detectors |33|, 34, |3^, |36| . 

At extremely-high energies (E > 5 x 10 17 eV), evidence was presented for neutral particles from 
the direction of Cygnus X-3 by the Fly's Eye |37j and by Akeno |3^[ groups, based on data taken 
during the periods 1981-1989 and 1984-1989, respectively. These data apparently indicated steady 
emission of neutral particles from Cygnus X-3 that was uncorrelated with the X-ray periodicity. 
The Haverah Park experiment, operating in the same energy range, and during much of the same 
period in time, found no evidence for such emission |3 
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The evidence from ground-based experiments for gamma-ray emission from Cygnus X-3 from 
1975 to 1990 is shown in Figure |l|. Here, the integral gamma-ray fluxes are plotted as a function 
of energy. Also shown is a single power law fit of the form E . The fact that the gamma-ray 
fluxes at widely varying energies could be approximately fit by a single power law was taken by 
some as evidence of a unified acceleration mechanism at the source. One should note, however, 
that all results shown in Figure |l] represent integral flux measurements by experiments incapable 
of accurately measuring differential fluxes. Since the detections were generally only marginally 
statistically significant, the reported fluxes equally represent the three to four standard deviation 
sensitivity of each instrument at a fixed energy. The fluxes therefore would naturally fall on 
an E _1 power law, if the sensitivities of the experiments scaled linearly with energy (which was 
approximately true for these first-generation experiments). It has also been pointed out that even 
if the source mechanism produced emission with a single power law form, the detected flux at 
Earth would have a significant dip between 10 3 and 10 4 TeV because of absorption of gamma-rays 
by the cosmic microwave radiation [f40| . 

Starting with the CYGNUS experiment in 1988 |41[| , a number of more sensitive ground-based 
experiments were unable to detect gamma-ray emission from Cygnus X-3, at levels significantly 
lower than the earlier reports. Upper limits on the flux were reported for experiments using both 



the atmospheric Cherenkov technique [42], as well as the ground- array technique [|4^, Ejj. Using 



parts of the eventual CASA-MIA detector, some of us reported limits for data taken in 1988-1989 



[45] and in 1989 pq] . The general trend of a "quiet" Cygnus X-3 continued into the early 1990's, 



although there were several publications claiming gamma-ray emission based largely on data that 



had been taken in the previous decade |47|, 48 



3.2 Hercules X-l 

Like Cygnus X-3, Hercules X-l is a powerful binary X-ray source [p3j. The X-ray emission is 
modulated on a time scale of 1.7 days which is assumed to result from the orbital motion of the 
binary pair. Unlike Cygnus X-3, Hercules X-l is not seen in radio, but has been observed for 



many years in the optical range [51| and a 5.8 kpc distance to the source has been determined 
[52]. In addition, Hercules X-l contains an X-ray pulsar with a period of 1.24 sec |j50| , but whose 
ephemeris is relatively poorly determined because of unpredictable variations in the spin-up rate 
Hercules X-l has not been detected by space-borne gamma-ray instruments. 
The first evidence from a ground-based observatory for gamma-ray emission from Hercules 



X-l was reported in 1984 by the Durham group using the atmospheric Cherenkov technique [54]. 
The reported gamma-ray emission came in the form of a short burst (~ 3 minute duration) that 
exhibited 1.24 sec periodicity. Following this report, additional pulsed emission was claimed by 



Cherenkov detectors operating at ultra-high energies (E > 500 TeV) [55] and at TeV energies 

mm. 



The most intriguing evidence for gamma-ray emission from Hercules X-l came from data 
taken in 1986 by three experiments. Data taken between April and July of 1986 by the Haleakala 



and Whipple |59j telescopes operating near 1 TeV, and by the CYGNUS experiment |60j in 
July of 1986 operating at energies above 50 TeV, all indicated evidence for gamma-ray emission 
from Hercules X-l in the form of short bursts of approximately 0.5 hr duration. In addition, 
the emission detected by each experiment exhibited a common periodicity near 1.2358 sec, which 
differed by a significant amount (~ 0.16% lower) from the known X-ray period. The data from the 
CYGNUS experiment was further puzzling because the events from the direction of Hercules X-l 
had a muon content that was similar to the cosmic ray background events, whereas gamma-ray 
events should have contained significantly fewer muons. Later, two groups with somewhat poorer 
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sensitivity presented additional evidence for gamma-ray emission from Hercules X-l at different 



times in 1986 at TeV || and 100 TeV II energies. 

Since the advent of upgraded and improved experiments in the early 1990's, the gamma-ray 
signals from Hercules X-l disappeared from the published literature. The Whipple group, using 
a more sensitive Cherenkov imaging technique, failed to detect emission from Hercules X-l, and 
found no statistically significant evidence for gamma-ray emission from Hercules X-l over a six 



year period, even including their data from 1986 [63]. An enlarged and improved CYGNUS 
experiment also failed to see gamma-rays from Hercules X-l in the period between 1987 and 1991 
Using data taken in 1989, we reported upper limits on the emission of gamma-rays from 



Hercules X-l using part of the eventual CASA-MIA experiment [46]. 
3.3 Theoretical Implications 

The many claims of very high energy gamma-ray emission from the binary systems Cygnus X- 
3 and Hercules X-l fueled great interest in the development of astrophysical models to explain 
such emission. There were also non-standard particle physics models put forward to explain the 
observations; these models will not be discussed here. 

For the case of Cygnus X-3, where the gamma-ray emission was generally observed with a 4.8 hr 
periodicity, the astrophysical models needed to incorporate the orbital dynamics of the binary 
system. Models in which an energetic pulsar alone served as the power source for the gamma-rays 



|]65| or in which accretion powered the gamma-rays [66] were proposed. These models generally had 
difficulty in producing gamma-rays at energies above 10 15 eV. More popular were a general class 
of models in which the gamma-rays were produced from the decays of 7r° mesons made in hadronic 



collisions [37, 68, 69]. The hadronic beam resulted from diffusive shock acceleration, perhaps near 
the neutron star, and possible beam targets included the atmosphere of the companion star or 
material in the accretion disk. Such "beam-dump" models were capable of explaining gamma-rays 
at ultra-high energies and were also able to accommodate the observed periodicities of the gamma- 
ray signals. Several authors recognized that in order to explain the ultra-high energy gamma-ray 
fluxes initially seen, the required luminosity of Cygnus X-3 would also be sufficient to account for 



a substantial fraction of the high energy cosmic ray flux [7C]. Hillas pointed out that if Cygnus 
X-3 consisted of a 10 17 eV accelerator with a luminosity of ~ 10 39 ergs/sec, only one such object 
like it would be required to explain the origin of cosmic rays above 10 16 eV f71Jl . 

Unlike Cygnus X-3, the gamma-ray emission from Hercules X-l was not seen to be correlated 
with the orbital motion of the binary system, but instead with the pulsar periodicity. This 
observation, along with evidence that the emission appeared in the form of short bursts, led 
naturally to models in which the pulsar itself was the power source. In such models, the gamma- 
rays were produced by the interaction of a charged particle beam with the accretion disk |f72| , 
73 [ . More difficult to explain were the 1986 observations of gamma-ray emission at a slightly 



shorter period than the X-ray period. The anomalous gamma-ray periodicity was explained by 
the presence of matter in the accretion disk which periodically obscured the gamma-ray interaction 



region J74J, |75j . 



In summary, although theoretical difficulties existed in explaining the apparent signals of 
gamma-rays from Cygnus X-3 and Hercules X-l, the signals were tantalizing because of the pos- 
sibility that they revealed important sources of the ultra-high energy cosmic ray flux. 



4 Experimental Procedure 
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Table 1: Size and makeup of CAS A- MIA experiment as a function of time for the data sample 
used in this analysis. A range of values indicates that the experiment was being enlarged during 
this period of time. Data taken after August 1995 are not used in this analysis. 





Epoch 




Enclosed Area (m 2 ) 


CASA Detectors 


MIA Counters 


Mar. 


1990 - Oct. 


1990 


108,900 


529 


512 


Oct. 


1990 - Apr. 


1991 


108,900-230,400 


529-1089 


512-1024 


Jan. 


1992 - Aug. 


1993 


230,400 


1089 


1024 


Aug. 


1993 - Aug. 


1995 


216,225 


1056 


1024 



4.1 The C ASA-MI A Experiment 

The CASA-MIA experiment is located in Dugway, Utah, USA (40.2° N, 112.8° W) at an altitude 
of 1450m above sea level (870g/cm 2 atmospheric depth). CASA-MIA consists of two major 
components: the Chicago Air Shower Array (CASA), a large surface array of scintillation detectors, 
and the Michigan Array (MIA), a buried array of scintillation counters sensitive to the muonic 
component of air showers. 

CASA consists of 1089 scintillation detectors placed on a 15 m square grid and enclosing an 
area of 230,400 m 2 . Construction on the array started in 1988 and a small portion (5%) of the 
experiment operated in 1989. A more substantial portion (~50%) of it was completed by early 
1990. Data collection with this portion started on March 1, 1990. Additional detectors were 
added in 1990 to complete the construction. 

MIA consists of 1024 scintillation counters located beneath CASA in 16 groupings (patches). 
The total active scintillator area is 2, 500 m 2 and the counters are buried beneath ~ 3.5 m of soil. 
This depth corresponds to a muon threshold energy of approximately 0.8 GeV. Parts of MIA were 
operational as early as 1987, 50% of the experiment was completed by early 1990, and the entire 
array was working by early 1991. The CASA-MIA experiment was turned off temporarily in 1991 
for repair due to lightning damage, but has operated essentially uninterrupted since that time. 
Table [l] summarizes the size and detector makeup of CASA-MIA as a function of time. 

Figure |2| shows a plan view of the experimental site. In addition to CASA-MIA, there are 
other installations at the same site. The other equipment used in this analysis is an array of 
five tracking Cherenkov telescopes. One telescope is located at the center of CASA-MIA, and 
the other four are 120 m away from the center along the major axes. Each telescope consists of 
a 35 cm diameter mirror which focuses Cherenkov radiation onto a single 5.1cm photomultiplier 
tube (PMT). The signals from the PMTs are digitized to record the amplitude and time of arrival 
of the Cherenkov wavefront at each telescope location. The shower direction is reconstructed by 
fitting the Cherenkov arrival times to a conical wavefront. 

A complete description of the CASA-MIA experiment can be found elsewhere 76]; here we 



briefly describe some aspects of the experiment that are relevant for this analysis. Each CASA 
station consists of four scintillation counters connected to a local electronics board. A station 
is alerted when at least two of the four counters fire within a 30 nsec window and it is triggered 
when at least three of the counters fire. If three or more stations trigger within a time period 
of approximately 3 //sec, an array trigger is said to have occurred. The array trigger rate de- 
pends on operating conditions (e.g. atmospheric pressure), but is ~20 Hz for the full CASA-MIA 
experiment. 

Upon an array trigger, the Universal Time (UT) is latched and recorded by either a GOES 
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Satellite Receiver Clock (1990-1993) with a precision of ±1 msec, or by a Global Positioning System 
(GPS) clock (1993-1995) with a precision of ±100nsec. For each array trigger, a command is 
broadcast to the array instructing each alerted CASA station to digitize and record its data, and 
a signal is generated to stop time-to-digital converters (TDCs) on each MIA counter. The TDCs 
have a range of 4 /xsec and a least bit precision of 4 nsec. The data from each CASA station consist 
of the arrival times and pulse-height amplitudes of the pulses from each scintillation counter, as 
well as the arrival times of pulses from the four nearest neighbor stations. The data from each 
MIA counter consist of the time of arrival of the array trigger relative to the passage of a muon 
through the counter. The CASA-MIA data and the Universal Time recorded as the result of an 
array trigger correspond to a single air shower event. 



4.2 Event Reconstruction 

We briefly summarize some of the important aspects of the CASA-MIA event reconstruction; full 



details can be found elsewhere [76]. The data from the experiment are accumulated in runs of six 
hours duration. All calibrations and offsets are determined for each run separately. At the start 
of a run, the timing constants associated with the CASA station electronics are calibrated by an 
internal oscillator. Timing constants are corrected for the effects of temperature by studying the 
constants over the span of a week. The CASA counter particle gains are determined for each run 
from the abundant cosmic ray air showers. The counter gains are found from the PMT amplitude 
distributions of those counters hit in stations with two out of four counters hit. A statistical 
correction of ~ 20% accounts for the fact that on average slightly more than one particle passes 
through a counter in this situation. The CASA cable and electronic delays are determined from 
the zenith angle distributions of the detected events. The relative delay between the CASA and 
MIA trigger systems is determined by centering the peak of the muon arrival time distribution 
relative to the position of the CASA trigger time. 

We estimate the shower core position by the location on the ground with the highest particle 
density. The total number of particles in the shower, or shower size, is determined by fitting the 
density samples obtained from the CASA stations to a lateral distribution of fixed form. The 
mean number of alerted CASA stations is 19 and the mean shower size is ~ 25,000 equivalent 
minimum ionizing particles. 

The shower direction is determined from the timing information recorded by CASA. The 
relative times between pairs of adjacent alerted CASA stations are determined. Each relative 
time gives a measure of the shower direction along one axis of the experiment. The times are 
weighted by an empirical function of the local particle density and distance to the shower core, 
and are fit to a wavefront which accounts for the conical shape of the shower front. The cone 
slope is approximately 0.07nsec/m. The shower direction in local coordinates is defined by two 
angles. The zenith angle, 9, is measured with respect to the vertical direction and the azimuthal 
angle, eft, is measured with respect to East in a counter-clockwise manner. 

In order to be confident of any astronomical results, it is essential to measure the angular 
resolution of the experiment. The resolution has two parts. The statistical part largely derives 
from the intrinsic fluctuations in the arrival times of the shower particles and from the timing 
resolution of the CASA counters and electronics. The systematic contribution derives primarily 
from the accuracies of the experiment survey and of the calculation of timing delays and offsets. 

The statistical contribution to the angular resolution is determined by three different tech- 
niques. First, on an event-by-event basis, we divide the array into two overlapping sub-arrays 
and compare the shower directions that are reconstructed by each sub-array. Using an air shower 
and detector simulation, we estimate the statistical correction required to derive the angular res- 
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olution from the sub-array direction comparison. Second, we compare the shower direction as 
determined by CASA with the direction determined by the five Cherenkov telescopes for events 
in which both CASA and the telescopes triggered. By statistically removing the angular resolu- 
tion of the telescopes from this comparison, we estimate the CASA resolution. Third, we have 



detected the shadow that the Moon casts in the cosmic rays [77|. For data taken between 1990 
and 1995, the Moon shadow is shown in Figure ||. By deconvolving the size of the Moon from 
the measured shadow, we obtain another estimate of the angular resolution. Figure || shows the 
resolution estimates from these different techniques. The agreement between the various methods 
is good, which allows us to determine a single parametric form for the resolution as a function of 
the number of alerted CASA stations. 

The systematic contribution to the angular resolution of the experiment has been checked 
by two different techniques. First, for data taken in coincidence with the tracking Cherenkov 
telescopes, we examine the angular difference between the directions determined by CASA and by 
the telescopes. The distribution of these differences indicates that the systematic offset between 
CASA and the telescope array is very small (< 0.1°). The alignment of the telescope array has 
been verified by the observation of a number of stars. A second check on the pointing accuracy of 
CASA comes from the Moon shadow. The center of the Moon shadow image is within 0.1° of the 
known position of the Moon. We conclude that the pointing uncertainty of CASA is negligible in 
comparison with the experiment's angular resolution. 

The muon content of the shower is determined from the data recorded by MIA. Since MIA 
records the times of counter hits over an interval of 4 [isec, it is sensitive to muons produced by 
showers arriving at any location of the array and from any direction. During the same time interval, 
MIA also records accidental counter hits produced by PMT noise and by natural radioactivity in 
the ground. The average number of accidental hits is approximately sixteen per event, while the 
average number of real muons associated with air showers is approximately nine per event. 

Real muons arrive within lOOnsec of the shower front arrival, while the accidental hits occur 
randomly over the 4 /isec interval. We greatly reduce the acceptance for accidental muon hits by 
narrowing the time window for accepting muons. The width of the window is determined from the 
distribution of muon times for each six hour run. We set the width to encompass 95% of the real 
shower muons; on average, it is ~ 150nsec. The position of the window is found on an event-by- 
event basis by means of a clustering algorithm. The algorithm searches for the cluster of three or 
more muons within an interval of 40nsec. In approximately 25% of the events, no cluster is found 
and the window position is placed at the center of the muon time distribution as determined for 
the entire run. As a result of tightening the time window for muon hit acceptance, the average 
number of real muons recorded is 8.5 per event, while the average number of accidentals is 0.63 
per event. 

The CASA-MIA data undergo several stages of processing and compression. In the most highly 
compressed format upon which this analysis is based, the data records are 26 Bytes per event and 
include the following information for each event: Universal Time (UT), number of alerted CASA 
stations, number of in-time muon hits, core location, arrival direction, shower size, and muon 
shower size (not used here). 

5 Analysis 
5.1 Data Sample 

The data used in this analysis were taken between March 4, 1990 and August 10, 1995, with a gap 
of 255 days in 1991. The experiment had usable data on 1627 days with the remainder of the days 
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Table 2: Quality cut efficiencies and event totals for the CASA-MIA data sample. The muon-data 
quality cuts are applied after the all-data quality cuts. The data sets (all and muon) are described 
in the text. 



Category 


All-Data Sample 


Muon-Data Sample 


Initial Event Total 


2087.8M 


1925. 8M 


Run Cut Efficiency 


0.935 


0.929 


Event Cut Efficiency 


0.986 


0.896 


Overall Efficiency 


0.922 


0.832 


Final Event Total 


1925. 8M 


1602. 7M 



lost largely because of power outages at the site and computer problems. The experiment has an 
instrumental deadtime of approximately 5.4% which is due to a number of effects, including data 
acquisition computer latency and the time needed to digitize the CASA station data. Calibration 
runs of approximate length of six minutes taken at the start of data runs, losses due to 8 mm tape 
failures, and downtime from array maintenance led to an additional reduction in the live time to 
a total of 1378.4 days (84.7% of the total). After the reduction and processing of the data, the 
final data sample consisted of 2.0878 x 10 9 reconstructed events. 

The size of the CASA-MIA data sample is unprecedented in air shower physics. To ensure data 
integrity, we impose a comprehensive set of data quality cuts. The cuts are tailored separately 
for the data sample in which we only use information from the surface array (all-data) and the 
sample in which we use information from both the surface and muon arrays (muon-data). For 
each of these samples, we make quality cuts on an event-by-event basis and on a run-by-run basis. 
Cuts are applied to runs and events only in the cases there there is evidence of an instrumental 
bias. The efficiencies of the cuts are summarized in Table ||. 

For the all-data sample, the run and event cuts have a combined efficiency of 92.2%, which 
yields a final sample of 1.9258 x 10 9 events. The most restrictive run cut requires a minimum 
fraction of the CASA stations to be working reliably and removes 2.2% of the data, largely because 
of instances in which isolated parts of the array failed. For the muon-data sample, the run cuts 
have an efficiency of 92.9%. A cut which requires a sufficient fraction of the muon counters to be 
working removes 4.8% of the data. The event cuts have an additional efficiency of 89.6%. The 
most restrictive event cut eliminates 3.2% of the events because they have no muon information 
due to deadtime of the MIA data acquisition system. The overall efficiency of the muon-data cuts 
is 83.2%, which yields a final sample of 1.6027 x 10 9 events. 

5.2 Gamma-Ray Selection 

From prior observations of Cygnus X-3 and Hercules X-l, we expect that gamma-ray fluxes, if 
present, will be small in comparison with the isotropic cosmic ray flux. Therefore, we need to 
enhance the presence of a possible gamma-ray signal by eliminating as many cosmic ray air showers 
as possible, while keeping a high fraction of the gamma-ray air showers. To do this, we select 
those showers with a reconstructed direction consistent with the position of the sources (within 
the angular resolution of the experiment) and with a muon content consistent with that expected 
from a gamma-ray primary. 
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Table 3: Angular search bin sizes and event fractions as a function of the number of CAS A alerts. 
The search bin is a circular region in equatorial coordinates whose radius is equal to 1.59 times 
the angular resolution. 



Alert Range 


Event Fraction 


Search Bin Radius 


3 - 10 


0.331 


2.45° 


11 - 15 


0.224 


1.88° 


16 - 20 


0.121 


1.40° 


21 - 30 


0.150 


1.05° 


31 - 40 


0.064 


0.78° 


41 - 60 


0.058 


0.60° 


> 60 


0.052 


0.41° 



5.2.1 Angular Search Bin 

We define a circular search bin whose size is based on the estimated angular resolution of the 
experiment. For a sufficiently large number of events, the bin which optimizes the signal-to- 
noise has a size equal to 1.59 times the angular resolution and contains 72% of the signal. The 
CASA-MIA angular resolution depends on the number of alerted CASA stations in an event, and 
therefore we use a variable-sized search bin which scales with the number of alerts. For simplicity, 
we use seven different bin sizes that range from 2.45° radius for showers with the least number of 
alerts, to 0.41° radius for showers with the largest number of alerts. These bin sizes are shown in 
Table ||, along with the fraction of events in each alert range. 

5.2.2 Muon Content 

Air showers created by gamma-ray primaries are expected to contain far fewer muons than showers 
initiated by cosmic ray nuclei. This expectation results because the cross section for photo- 



pion production is much smaller than the cross section for electron-positron pair production [78]. 
Therefore, the interaction of a high energy gamma-ray in the atmosphere is much more likely to 
produce an electromagnetic cascade in the atmosphere than it is to create a hadronic cascade. 
Conversely, cosmic ray nuclei preferentially interact to create hadronic cascades. Showers initiated 
by gamma-rays are thus expected to contain far fewer hadrons than those initiated by cosmic rays. 
Since air shower muons are predominantly produced from the decays of pions and kaons in the 
hadronic cascade, gamma-ray air showers should contain far fewer muons as well. Simulations have 
been done to estimate the muon content of air showers @ O ur own simulation indicates 
that an air shower initiated by a 100 TeV gamma-ray contains, on average, 3-4% of the number 
of muons in a shower initiated by a proton of the same energy. 

The muon content of showers should in principle be a powerful tool in rejecting cosmic ray 
background events. In our experiment, the rejection capability is limited by the collection area 
of the muon array and, to a lesser extent, by the presence of a small amount of accidental muon 
hits. The muon array (MIA) is significantly larger than any other air shower muon detector built 
to date, but its active area still corresponds to only ~ 1% of the enclosed area of the experiment. 
As shown in Figure |, the average of the distribution of the number of in-time muons is ~8.5, but 
the shape of the distribution is such that its mode is three, and a substantial fraction of events 
have zero muons. In Figure || we also show the estimated number of muons for showers initiated 
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Table 4: Quantities associated with the selection of muon-poor events. Muon-poor events are 
those having a relative muon content, r^ (defined in the text), less than a cut value. The cut 
values are given in the second column, and the third and fourth columns give the efficiencies for 
passing the cut for gamma-ray signal events and for hadronic background events, respectively. 
The fifth column gives the quality factor, Q, or the improvement in flux sensitivity from making 
the cut. 



Data Set 


r^ cut Value 


Signal e 


Background e 


Q 


All 


-0.75 


0.72 


0.0600 


2.94 


< 10 Alerts 


-0.50 


0.69 


0.1644 


1.70 


> 10 Alerts 


-0.75 


0.76 


0.0362 


3.99 


> 40 Alerts 


-1.00 


0.71 


1.77 x 10~ 3 


16.9 


> 80 Alerts 


-1.00 


0.77 


0.67 x 10" 3 


29.7 



by gamma-rays, including the contribution from accidental muon hits. For gamma-ray showers, 
we expect, on average, 0.28 real muons per event and 0.63 accidental muons per event. 

In order to enhance a possible gamma-ray signal, we wish to select muon-poor events, i.e. events 
that have fewer muons than the average expected number. To do this, we make the assumption 
that any gamma-ray signal in the data is much smaller than the flux of cosmic rays. We can 
therefore use the muon information from the detected events to describe the muon content of the 
background, and our simulation to describe the muon content of the gamma-ray signal. 

The number of muons in a shower depends on a number of observable quantities, for example, 
the number of alerted CASA stations, shower zenith angle, and core position. We develop a 
parameterization for the average number of muons as a function of these quantities by examining 
a large ensemble of actual showers. We then determine the relative muon content of a specific 
shower by comparing the observed muon number, (n jU ) b s , to the expected number of muons, 
< n M > 0X p, for showers having similar zenith angles, core positions, and numbers of alerts. The 
relative muon content, r M , is defined by: 



= Log 



10 



(n 



fi Jobs 



(1) 



Figure ^ shows the distributions of r^ for observed events and for simulated gamma-ray events. 
Muon-poor events are defined as those having r^ values less than some cut value. The position of 
the cut is chosen to reject as many background events as possible, while keeping a high fraction 
of the gamma-ray events. The cut value depends weakly on the number of CASA alerts because 
the separation between the signal and background r^ distributions improves as the showers get 
larger. 

Table || shows the r^ cut values for various samples of data along with the fractions of signal 
and background events retained, and the sensitivity improvement achieved from making a cut. 
For the entire data set, the sensitivity is improved by a factor of 2.94 by cutting on the shower 
muon content. The quality factor increases to 29.7 for events having more than 80 alerted CASA 
stations. 
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5.3 Background Estimation 



We select gamma-ray candidate events {on-source events) based on their reconstructed arrival 
direction in equatorial coordinates (right ascension, a, and declination, 5) and on their muon 
content. In order to derive the significance of a possible gamma-ray signal, we need to determine 
the expected number of background cosmic ray events (off-source events) that would arrive from 
the same direction in the sky as the source and would have a similar muon content as gamma-ray 
events. Again, we make the assumption that the detected air showers are predominantly caused 
by background cosmic ray events. We thus use the detected events themselves to estimate the 
expected background. 

A common method to estimate the expected number of background events is to use off-source 
bins having the same declination as the source, but having different right ascension values. This 
method, which assumes a uniform experiment exposure over declination, was satisfactory for 
earlier smaller experiments. However, given our large event sample, this technique is not reliable 
for CASA-MIA because of small, but non-negligible, systematic biases (e.g. diurnal variations). 
For the CASA-MIA data sample, a source at a declination of 40° occupies an angular bin with 
~ 1.8 x 10 6 events. The fractional statistical uncertainty corresponding to one standard deviation 
in the number of events is 0.075%. In order to accurately estimate the number of background 
events, the relative systematic uncertainty must be well below this level. As a result, an accurate 
and robust way to determine the expected background is needed. Several methods have been 
developed by other groups 37, 81] and the method that we use is similar to these. 

The detection rate of an air shower array triggering on cosmic rays is determined by the 
properties of the cosmic ray flux and by the properties of the array itself. Assuming that the cosmic 
ray parameters do not change with time, any variation in the detection rate is caused only by 
changes in the detector or in the atmospheric conditions. Over short intervals of time (~ lhour), 
the relative detection efficiency as a function of the shower direction in local coordinates, (0,(f>), 
is largely determined by the array geometry (placement of detectors, uniformity of terrain, etc.) 
and is almost constant, and the time variation of the detector response may be estimated from the 
trigger rate. Therefore, we separate the detection rate per unit solid angle in local coordinates, 
N(9, (j), t), into two terms: 



N{6,(/>,t) = D{6A)-R{t) , (2) 

where D(9, <j>) is the efficiency per unit solid angle of detecting a shower from a given direction in 
the sky, and R(t) is the trigger rate as a function of time. The factor D(0, 4>) is determined by 
maps made from the arrival directions of cosmic ray showers over given periods of time. The time 
dependent term, R(t), is determined from the arrival times of the actual events. The expected 
number of events for a given bin in the sky is then determined by integrating N(6, (ft, t) over the 
time interval in question. To determine the expected number of events for a bin in equatorial 
coordinates, (a, 5), we integrate N(0,cp,t) over the time interval and over local coordinate space. 

More explicitly, the background estimation is done by the following procedure. For intervals of 
4, 200 sec, we accumulate the arrival directions of cosmic ray events into 2,700 bins segmented in 
local coordinate space (30 bins in 8, 90 bins in <p). We use the binned data to construct maps of the 
relative acceptance of any point in the sky over this time interval. Separate maps are calculated 
for each data sample used in the source search (e.g. all-data and muon-poor data). To generate 
simulated background data, we discard the directional information of an event and associate the 
event time with a local coordinate direction obtained by sampling from the appropriate sky map. 
We then compute artificial values for the equatorial coordinates and determine if this simulated 
event falls into a search bin of a source. By sampling more than once from the sky map for each 
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event time, we increase the statistics on the simulated data sample. Negligible systematic bias 
is introduced by such oversampling. For this work, we oversample by a factor of ten, an amount 
that is limited only by computational resources. 

We have checked that our background estimation method is free from bias by comparing the 
detected numbers of events in an angular bin to our expected number for bins that do not contain 
Cygnus X-3 and Hercules X-l. For each bin, we compute the statistical significance of any excess 
or deficit in the number of detected events relative to the number we predict. The distribution 
of these significances is in close agreement with that expected from statistics, which, because of 
background oversampling, is dominated by the statistical uncertainty on the number of detected 
events. 



5.4 Energy Response 

Air shower arrays trigger on the shower size, i.e. the number of charged particles in the shower at 
ground level. For each shower, we determine a shower size from the particle densities measured 
in the CASA stations. For astrophysical interpretations of flux measurements or flux limits, 
however, it is necessary to translate from the measured shower parameters (size and zenith angle) 
to an estimate of the energy of the primary particle. Since there are large fluctuations in shower 
size for showers initiated by particles at fixed energy, it is difficult for air shower experiments 
to measure accurately differential primary spectra. Traditionally, therefore, flux measurements 
have been quoted as integral intensities above a fixed energy point. Although to some degree the 
energy value at which to quote the intensity is arbitrary, we desire to use an energy at which 
the experiment has a significant degree of sensitivity. We chose to quote flux measurements at 
the median energy of the experiment which reduces the dependence of the flux on the assumed 



spectral index |45|, |85| . 



We estimate the energy response of the experiment by the constant intensity method, which 
has been used by other experiments as well as by our own group [^ ]. The constant intensity 
procedure is described in more detail elsewhere [84]. Briefly, we determine the relationship between 
shower size and energy by comparing the detected flux of showers above a given size to an assumed 
form of the all-particle cosmic ray spectrum. The comparison is done on a run-by-run basis to 
account for changes in the detector response. The cosmic ray flux is derived from measurements 



made by other space-borne [gg, 87] and ground-based Q] experiments. The assumed integral 
cosmic ray intensity above 100 TeV is 6.57 x 10 -9 particles cm~ 2 s _1 sr _1 . 

We use the relationship between energy and shower size to determine the most likely energy 
for each shower coming from the direction of Cygnus X-3 or Hercules X-l in a angular bin of fixed 
radius. The medians of the energy distributions determine the median energies for cosmic ray 
particles from the direction of Cygnus X-3 and Hercules X-l that would trigger the experiment 
and pass all selection criteria. By normalizing our energy scale to the cosmic ray flux, we make 
the assumption that the primary particle has the same spectral index as the detected cosmic rays. 
This assumption is reasonable when dealing with sources like Cygnus X-3 and Hercules X-l in 
which there are are no well established measurements of spectral indices. The median energy of 
particles from the direction of Cygnus X-3 is 114 TeV, and from Hercules X-l, it is 116 TeV. Since 
the difference in the energies for the two sources is negligible, we report our measurements at a 
common energy of 115 TeV. 



5.5 Search Strategy 

We carry out searches for particle emission from a particular source by comparing the number of 
events found within a circular angular bin around the source to the number of events estimated by 
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Table 5: Data samples selected by integral cuts on the number of CASA alerts. The cut values 
are given in the first column and the fractions of events surviving the cut (and within the angular 
search region) are shown in the second column. The median energies for events coming from either 
Cygnus X-3 or Hercules X-l are listed in the third column. 



Alert Cut 


Event Fraction 


Median Energy 


None 


100.00% 


115 TeV 


< 10 


62.87% 


85TeV 


> 40 


0.58% 


530 TeV 


> 80 


0.09% 


11 75 TeV 



our background procedure. The angular bin sizes vary as a function of the number of alerted CASA 
stations, as itemized in Table ||. Source positions (J1992) are taken to be (a, 5) = (308.04°, 40.93°) 
for Cygnus X-3, and (a, 5) = (254.39°, 35.35°) for Hercules X-l. 

Separate searches are made based on particle type and energy. By using the all-data sample, 
we are sensitive to any type of neutral particle that would create air showers. With the muon-poor 
sample, we are specifically sensitive to the emission of gamma-rays. We carry out three separate 
searches with various integral cuts on the number of alerted CASA stations, in addition to a search 
with no cuts. This procedures takes advantage of the correlation between primary energy and size 
(as represented by the number of alerts), and improves our sensitivity to possible emission that 
might be present at either low or high energies. The data samples selected by cutting on the alert 
number and their corresponding median energies are shown in Table ||. 



6 Results 

We search for evidence of neutral (gamma-ray or other) particle emission from Cygnus X-3 and 
Hercules X-l. Separate searches are carried out for steady and transient emission from either 
source. In addition, we search for periodic emission from Cygnus X-3 at the 4.8 hr X-ray periodicity 
and for emission from Cygnus X-3 that was coincident with the occurrence of large radio flares. 
No compelling evidence for emission from either source is found for all the different searches, and 
consequently we set upper limits on the fluxes of particles from the sources. 



6.1 Steady Emission 

The numbers of on-source and background events for the various searches from Cygnus X-3 are 
shown in Table [6|. The results from similar searches carried out on Hercules X-l are shown in 
Table [7[ For each search, we also calculate the statistical significance of any excess or deficit 
in the number of events observed relative to background by the prescription of Li and Ma [88], 
using an oversampling factor of 10. No significant excess is observed for any search from either 
source. Therefore, for each search, we calculate an upper limit, Nqq, on the number of excess 
events from the source at the 90% confidence level [^9], |9(| . Each N$q value is converted to a limit 
on the fractional excess of events from the source, /go, by dividing by the estimated number of 
background events, which is assumed to represent the background cosmic-ray level. Since the /go 
values are independent of the absolute flux normalization, they are useful in comparing results 
between different experiments. 
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Table 6: Steady emission search results for Cygmis X-3 using the all-data sample (top) and muon- 
poor sample (bottom). The number of events observed on-source and the number expected from 
background are given in the second the third columns, respectively. The fourth column gives the 
statistical significance of any excess or deficit. The 90% c.l. upper limit on the number of excess 
events, Ngo, and the upper limit on the fractional excess, /90, are given in the last two columns. 
The methods used to calculate statistical significances and upper limits are outlined in the text. 
The data samples at 85TeV, 530 TeV, and 1175 TeV are subsets of the data sample at 115 TeV. 



All-Data Sample 



Energy 


On-Source 


Background 


Signif. 


^90 


/90 




85 TeV 


1119469 


1119987 


-0.48(7 


1502.1 


1.34 x 10" 


-3 


115 TeV 


1780594 


1781479 


-O.660- 


1774.9 


9.96 x 10" 


-4 


530 TeV 


10286 


10235 


+0.49(7 


205.3 


2.01 x 10" 


-2 


11 75 TeV 


1583 


1580 


+0.08(7 


68.9 


4.36 x 10" 


-2 


Muon-Poor Sample 


Energy 


On-Source 


Background 


Signif. 




/90 




85 TeV 


149676 


149863 


-0.57(7 


548.1 


5.90 x 10" 


-4 


115 TeV 


121409 


121594 


-0.37(7 


485.4 


3.28 x 10" 


-4 


530 TeV 


20 


21.0 


-0.21(7 


8.2 


9.47 x 10" 


-4 


11 75 TeV 


1 


0.6 


+0.44a 


3.5 


2.67 x 10" 


-3 



Table 7: Steady emission search results for Hercules X-l using the all-data sample (top) and 
muon-poor sample (bottom). 



All-Data Sample 



Energy 


On-Source 


Background 


Signif. 


^90 


/90 




85 TeV 


1058904 


1057583 


+ 1.12(7 


2738.1 


2.59 x 10" 


-3 


115 TeV 


1681708 


1681392 


+0.23(7 


2387.6 


1.42 x 10" 


-3 


530 TeV 


9579 


9532 


+0.46(7 


196.5 


2.06 x 10" 


-2 


11 75 TeV 


1419 


1459 


-0.98(7 


44.0 


3.02 x 10" 


-2 


Muon-Poor Sample 


Energy 


On-Source 


Background 


Signif. 


^90 


/90 




85 TeV 


139580 


139670 


-0.24(7 


577.1 


6.57 x 10" 


-4 


115 TeV 


113360 


113244 


+0.37(7 


643.4 


4.62 x 10" 


-4 


530 TeV 


14 


16.8 


-0.67(7 


6.3 


7.96 x 10" 


-4 


11 75 TeV 





0.5 


-0.98(7 


2.3 


1.90 x 10" 


-3 
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Table 8: Flux limits from searches for steady emission from Cygnus X-3 (top) and Hercules X-l 
(bottom). The second and third columns give the 90% c.l. upper limit on the integral flux of any 
neutral or gamma-ray particles from the source, respectively. The units of flux are particles cm -2 
sec -1 . 



Cygnus X-3 


Energy 


<Z> N (E) 


$7 ( E ) 




115 TeV 


2.20 x 10" 14 


6.26 x 10- 


-1S 


530 TeV 


1.43 x 10~ 15 


1.21 x 10 - 


-16 


1175 TeV 


1.04 x 10" 15 


5.19 x 10- 


-17 


Hercules X-l 


Energy 


$N (E) 


<*>7 (E) 




115 TeV 


3.04 x 10" 14 


8.55 x 10" 


-IS 


530 TeV 


2.87 x 10" 15 


9.75 x 10- 


-17 


1175 TeV 


6.91 x 10" 16 


3.56 x 10" 


-17 



Figure [?] shows scans in right ascension for bands of declination centered on Cygnus X-3 for the 
all-data and muon-poor samples. No significant excess above background is seen in either sample 
for the bin containing Cygnus X-3. The background estimation agrees well with the data in the 
off-source region. Similar scans for Hercules X-l are shown in Figure)^, and again the background 
estimation agrees well with the observed data and no excesses are seen. 

6.1.1 Flux Limit Calculation 

In the absence of a statistically significant excess from either Cygnus X-3 or Hercules X-l, we 
set upper limits on the flux of particles from each source. Separate limits are set for neutral and 
gamma-ray primaries. For gamma-ray primaries, the 90% c.l. upper limit, <J> 7 (-E), on the integral 
flux is calculated from the measured fractional excess by normalizing to the cosmic ray flux: 

M*0 = ^J(E). (3) 

Here, Cl is the mean solid angle used in the search, e is the fraction of events that would pass cuts 
and fall into the search bin, J(E) is the integral cosmic ray intensity above energy E, and i? 7 is 
a factor which accounts for the relative trigger efficiency for gamma-rays as opposed to cosmic 
rays. The value of Q ranges from 5.74 x 10 _3 sr for the lowest energy data set to 1.60 x 10 _4 sr 
for the highest energy data set. The e factor accounts for the fraction of gamma-rays that would 
end up in the angular search bin (0.72) and the fraction that would pass the muon-poor selection 
criterion (Table ||). The value of was determined by Monte Carlo simulations to be 1.6. 

To determine an upper limit on the integral flux of any neutral particle from a source, <&n(E), 
we use Eq. ||, except e is now 0.72 and i? 7 is 1.0. In this calculation, we assume that the neutral 
particle would interact in the atmosphere to create air showers in a similar manner to cosmic rays. 

Table |8] gives the flux limits obtained from the various searches for steady emission from the 
two sources. Limits are not calculated for the data samples with median energies of 85 TeV because 
these are not integral energy samples. 
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6.2 Transient Emission 



We search for transient emission of particles from Cygnus X-3 and Hercules X-l on daily (single 
transit) time scales. For each transit of the source, we compare the number of on-source events to 
the number of expected background events and calculate a significance based on the prescription 
of Li and Ma [pif ]. We require the live time fraction during the transit to be at least 0.20 to 
remove transits in which the experiment was operational for only a small fraction of the time. 
For each source, we make separate studies of the transit significances for the all-data and muon- 
poor samples, corresponding to possible emission from any neutral and gamma-ray particles, 
respectively. 

For Cygnus X-3, the number of good transits in the all-data sample is 1500. In the muon-poor 
sample, it is 1291. The distributions of significances for the two samples of Cygnus X-3 transits 
are shown in Figure |9[ Each distribution agrees well with a Gaussian distribution of mean zero 
and unit width. There is no evidence for any excess of events at high values of significance (either 
positive or negative). 

For Hercules X-l, there are 1492 good transits in the all-data sample, and 1271 good transits 
in the muon-poor sample. The significance distributions for Hercules X-l are shown in Figure [l(| 
and again, no evidence for significant excesses exists. 

Based on the lack of statistically significant excesses, we place limits on the daily fluxes of 
neutral and gamma-ray particles from Cygnus X-3 and Hercules X-l. These limits are calculated 
by a similar procedure as used for the steady searches. The limit values depend on the actual 
statistical significance of the search on a given day, and also on the epoch of data taking. As 
shown in Table [l], the size of the experiment has changed with time, and the sensitivity changed 
accordingly. In Table ||, we give typical daily flux limits for the two sources for different epochs of 
the experiment. Since the numbers of events detected per transit are the same for the two sources 
to within 5%, the limits for Cygnus X-3 and Hercules X-l are virtually identical. Typical 90% c.l. 
limits on the integral flux using the full experiment are &n(E > 115 TeV) < 9.7 x 10~ 13 neutral 
particles cm -2 sec" - 1 and <J> 7 (£ > 115 TeV) < 2.0 x 10~ 13 photons cm 2 sec . 

We have also carried out searches for transient emission on the shorter time scale of 0.5 hr. 
Here, we compare the number of events observed on-source to the expected background level for 
ten 0.5 hr time intervals on either side of the time of source culmination. The typical number of 
on-source events, and therefore the flux sensitivity, depends strongly on the source zenith angle. 
For example, for an overhead source near culmination, the experiment observes ~ 175 events per 
0.5 hr, whereas at four hours from culmination the rate is ~ 15 events per 0.5 hr. Regardless of 
the rate, for each 0.5 hr interval, we calculate the significance in the number of on-source events 
relative to the background and combine all such significances into a single distribution. The 
resulting significance distributions are consistent with those expected from background processes 
for both sources in both the all-data and muon-poor samples. The typical 90% c.l. upper limits 
on the fluxes from either source are &n(E > 115 TeV) < 3.1 x 10~ 12 neutral particles cm -2 sec -1 
and <S>^(E > 115 TeV) < 7.1 x 10~ 13 photons cm 2 sec 1 for 0.5 hr periods within one hour of 
culmination. 



6.2.1 Cygnus X-3 Radio Flares 

We study showers from the direction of Cygnus X-3 during the occurrence of large radio flares at 
the source. We define large flares as those times when the radio output at 8.3 GHz exceeded 2 
Jy, a level which is two orders of magnitude above the typical quiescent level. During the period 
of CASA-MIA operations, there were six large flares, as listed in Table |l2|, 13|. 
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Table 9: Typical daily upper flux limits (90% c.l.) for emission of neutral and gamma-ray parti- 
cles from Cygnus X-3 and Hercules X-l. The flux limits are calculated for two different epochs 
assuming the same number of on-source events as off-source. The third column gives the typical 
number of events observed on-source during the different epochs. Epoch I corresponds to March 
1990 to October 1990. Epoch II corresponds to January 1992 to August 1993. Flux limits for the 
remaining periods of time of operation are close to those for Epoch II. Units of flux are particles 
cm~ 2 sec -1 . 



Epoch I 



Energy 


Particle 


Events 


/go 


^ daily (E) 


115 TeV 


Any 


575 


0.071 


1.6 x W' 1 ' 2 


530 TeV 


Any 


3.2 


1.34 


1.9 x 10~ 13 


1175 TeV 


Any 


0.56 


4.11 


9.6 x 10~ 14 


115 TeV 


7-ray 


45 


0.026 


3.6 x 10~ 13 


Epoch II 


Energy 


Particle 


Events 


/go 


^ daily (E) 


115 TeV 


Any 


1450 


0.044 


9.7 x 10" 13 


530 TeV 


Any 


8.1 


0.76 


1.1 x 10~ 13 


1175 TeV 


Any 


1.3 


2.56 


5.9 x 10~ 14 


115 TeV 


7-ray 


96 


0.015 


2.0 x 10~ 13 



Table 10: Large radio flares of Cygnus X-3 from 1990 to 1995, coincident with the operational 
time of CASA-MIA. The flare number is an arbitrary index used for this work. The peak radio 
flux values (8.3 GHz) come from 12, O]. The March 1994 flare was actually a prolonged event 
that extended for the ten days following March 1, 1994. 



Flare 


Date 


Peak Flux (Jy) 


1 


Aug. 15, 1990 


7.5 


2 


Oct. 05, 1990 


10.2 


3 


Jan. 21, 1991 


14.8 


4 


Sep. 04, 1992 


4.1 


5 


Feb. 20, 1994 


4.9 


6 


Mar. 09, 1994 


5.2 
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Table 11: CAS A- MI A search results for emission from Cygnus X-3 near the time of large radio 
flares. The flare numbers are defined in Table [l(| The — and + designations refer to the days 
preceeding and following the flare day, respectively. The significances (columns 4 and 8) are 
standard deviation values calculated using the prescription of Li and Ma fl88|| . The last two 
columns give the 90% c.l. upper limits on the integral flux above 115 TeV of any neutral particle 
and gamma-rays, respectively, in units of 10 -12 particles cm -2 sec -1 . Entries having only a dash 
indicate the absence of any usable data. 
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We examine the daily significances for Cygnus X-3 on the day of each large flare, as well as 



on the day preceeding and following each flare. Table 11 lists the numbers of observed events, the 
expected background, and the Li-Ma significances for the examined days. There is no compelling 
evidence for any statistical excess in the observed number of events from Cygnus X-3 for either 
the all-data or muon-poor samples. On one day (Feb. 20, 1994) the Li-Ma significance is 2.26cr 
for the all-data sample. The probability that we would get a day with this level of significance or 
greater is 28.2% after accounting for the fifteen days in which we searched. In addition, on this 
same day, there is no evidence for any excess in the muon-poor data, while we would expect the 
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statistical significance to increase by a factor of 3.3 if it were due to a gamma-ray signal. Table 
also lists the derived fractional excess values, /go, as well as the upper limits to the integral flux 
of neutral or gamma-ray particles from Cygnus X-3 during the flares. 

6.3 Periodic Emission 

Several previous observations of Cygnus X-3 claimed evidence for steady emission correlated with 
the 4.8 hr X-ray periodicity of the source. For this reason, we carry out a search for such emission 
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Table 12: CASA-MIA search results for 4.8 hr periodic emission from Cygnus X-3. Flux limits are 
given for selected phase intervals in which earlier experiments had reported detections. Columns 
3 and 4 give the 90% c.l. upper limits to the integral flux of neutral and gamma-ray particles, 
respectively, in units of particles cm -2 sec -1 . A blank entry corresponds to a data set having 
insufficient data with which to calculate a limit. 



Phase Interval 


Energy 






d> 7 (£) 




0.2 - 0.3 


115 TeV 


8.9 x 10- 


14 


2.3 x 10- 


14 




530 TeV 


4.5 x lO - 


15 


6.9 x 10- 


16 




1175 TeV 


3.5 x 10" 


15 






0.6 - 0.7 


115 TeV 


1.4 x 10- 


13 


3.5 x 10" 


14 




530 TeV 


3.8 x 10" 


15 


3.4 x 10- 


16 




1175 TeV 


6.5 x 10" 


15 







using the entire CASA-MIA data set. The event arrival times (UT) are corrected to the barycenter 
of the solar system using the JPL DE200 planetary ephemeris |)l| . The corrected times are folded 
with the 4.8 hr X-ray ephemeris of van der Klis and Bonnet-Bidaud |pj. A slight correction is 
made for newer X-ray data from the ASCA satellite, as reported by Kitamoto et al. Q. Each 
event is then assigned a phase value in the interval (0,1) representing the fraction of a period that 
the event is from the X-ray minimum. The phase values are accumulated in twenty bins of 0.05 
phase units each for both the on-source and generated background events. 

Figure 11 shows the 4.8 hr periodicity distribution of events from the direction of Cygnus X-3 
for the all-data and muon-poor samples. Also shown is the phase distribution expected from the 
background events. No compelling excesses are seen at any particular phase interval for either 
sample. We carry out similar periodicity analyses using data at higher energies selected by the 
number of alerted CASA stations. These searches also do not indicate any significant excesses at 
any phase interval. In Table 12, we list flux limits for the various searches at the phase intervals 
(0.2,0.3) and (0.6,0.7). These intervals were ones in which numerous earlier experiments had 
reported detections. 



7 Comparison with Other Results 

As described earlier, the many detections of Cygnus X-3 and Hercules X-l by experiments oper- 
ating between 1975 and 1990 varied greatly in their characteristics. Some results were steady and 
some were episodic; some exhibited apparent periodicity and others did not. We thus choose to 
compare the results of this work to a generalized picture of the earlier results and to more recent 
work. 



7.1 Cygnus X-3 

In Figure [l^, we plot the flux limits reported here on steady emission from Cygnus X-3. We 
also show published results from other experiments using data taken at times which overlap our 
observation period. The other results come from the Tibet air shower array in Yangbajing, China 
p|, the CYGNUS array in New Mexico, USA §1, the EAS-TOP array at Gran Sasso, Italy 



and the HEGRA experiment on the Canary Island La Palma 95]. We do not show earlier results 
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from data taken by a portion of our experiment in 1989 46] or the results of our all-sky survey 



for northern hemisphere point sources using data taken in 1990-1991 [S3]. In Figure |l3|, we show 
a similar comparison of the limits on the fractional excess of events from Cygnus X-3 relative to 
the cosmic ray background. 

The data from the recent experiments are consistent; no steady emission of ultra-high energy 
particles (gamma-ray or otherwise) has been detected from Cygnus X-3 at levels which are con- 
siderably lower than earlier reports. At TeV energies, the results from the Whipple Telescope 
| |96| | are also considerably lower than the earlier reports. The limits presented here are a factor 
of 130 lower at 115 TeV, and a factor of 900 at 1175 TeV, than the spectrum plotted in Figure |]. 
Our results are also inconsistent with emission reported by a smaller experiment using data taken 
during a time that overlapped our observations |p7fl . 

The limits presented here on transient emission from Cygnus X-3 are lower than, but in 
agreement with, those reported by other air shower experiments. There have been no compelling 
reports of transient emission of gamma-rays from Cygnus X-3 over the period 1990-1995, including 
during large radio flares from the source. There was an observation of underground muons from 



the direction of Cygnus X-3 during the January 1991 radio flare [37|. The reported flux for this 
observation was 7.5 x 10~ 10 muons cm -2 sec -1 , for muon energies above 0.7 TeV. If the muons 
were produced in air showers by the interaction of a hypothetical neutral particle from Cygnus 
X-3, we would expect a typical neutral particle energy of ~ 10 TeV |98[j . Assuming that the 
particle spectrum continues to energies detectable by CASA-MIA (and conservatively using a soft 
spectrum comparable to the cosmic rays), one derives an expected flux of ~ 10 -11 particles cm -2 
sec -1 for energies above 115 TeV. This flux is a factor of 5 to 10 above the flux limits set by 



CASA-MIA on the emission of any neutral particle during the January 1991 flare (Table 11). 

We have also shown that there is no evidence for 4.8 periodic emission from Cygnus X-3. This 
result is consistent with reports by other experiments over the same period of time. The limits on 
pulsed gamma-ray emission presented here for the phase intervals of 0.2-0.3 and 0.6-0.7 (Table |l2|) 
are lower at 115 TeV, and considerably lower at 530 TeV, than the fluxes predicted by a recent 
theoretical paper pS| ]. 

7.2 Hercules X-l 

The limits on steady emission of gamma-rays from Hercules X-l presented here are in agreement 



with those from other experiments, as shown in Figure 14. Gamma-ray emission from Hercules 
X-l was typically seen by earlier experiments as transient emission over short time scales (e.g. 
the 1986 outbursts). We have no evidence for such emission over the entire period 1990-1995. 
In Figure |15L we compare the daily event totals observed by CASA-MIA from the direction of 



Hercules X-l to the total expected assuming the flux of an earlier reported outburst [60]. Clearly, 
no evidence for emission at even much weaker levels than this outburst is seen during this time. 
The flux reported in Ref. [6C] was ~ 2 x 10 -11 particles cm~ 2 sec -1 for minimum energies of 



100 TeV. This flux is about a factor of 45 larger than the typical limits placed by CASA-MIA 
during the early part of operations and about a factor of 80 larger than the typical daily gamma- 
ray limits placed by the full CASA-MIA experiment (Table |9|). Since we have no evidence for 
transient emission from Hercules X-l, we choose not to carry out a periodicity analysis based on 
the X-ray pulsar period of 1.24 sec. 
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8 Conclusions 



We have carried out a high statistics search for ultra-high energy neutral and gamma-ray particle 
emission from Cygnus X-3 and Hercules X-l between 1990 and 1995. We have no evidence for 
steady or transient emission from either source, and for Cygnus X-3, we have no evidence for 4.8 hr 
periodic emission or emission correlated with large radio flares. These results are in agreement 
with those from other experiments operating during the same period of time, but are in stark 
contrast to earlier (1975-1990) reports. 

The apparent disappearance of Cygnus X-3 and Hercules X-l from the ultra-high energy 
gamma-ray sky can be interpreted in two ways. An optimistic view Q is that the earlier re- 
sults indicated the presence of ultra-high energy gamma-rays (or particles) from Cygnus X-3 and 
Hercules X-l, and that the sources, which are episodic on long times scales, are now dormant. 
A more pessimistic view is that the earlier reported detections were largely, if not entirely, sta- 
tistical fluctuations, and that no compelling evidence exists for ultra-high energy gamma-rays 
from any astrophysical source. We point out that an earlier all-sky survey using a portion of our 
data sample indicates that the northern hemisphere does not contain any steady point sources 
of gamma-rays with fluxes comparable to those reported from X-ray binaries in the 1980's 



We have presented an update on this analysis at a conference | 100 | which are consistent with the 
absence of bright 100 TeV gamma-ray point sources. We are in the process of completing a final 
all-sky survey on the five year C ASA-MI A data sample. 

The pessimistic interpretation of the ultra-high energy point source question, if correct, high- 
lights the difficulties in detecting gamma-rays from sources at other (high) energies and in detecting 
neutrinos as well. In addition, without compelling evidence for high energy particle acceleration 
at point sources, the difficulties in explaining the origins of cosmic rays above 10 14 eV remain. 
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Figure 1: Published results from ground-based experiments indicating evidence for gamma- 
ray emission from Cygnus X-3 during the period 1975-1990. The circles indicate results from 
atmospheric Cherenkov telescopes, the squares show data taken by air shower arrays, and the 
triangle indicates the result from the Fly's Eye experiment. The dashed curve is an approximate 
power law fit to the data with a slope of —1.1. Not shown in this figure are several upper limits 
to the flux of gamma-rays from Cygnus X-3 during this same epoch. The two points at extremely 
high energy (5 x 10 5 TeV) have been slightly displaced from each other for clarity. 
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Figure 2: Plan view of the CASA-MIA experiment (Dugway, Utah, USA). Small squares indicate 
the 1089 scintillation detectors of the CASA surface array. Sixteen large rectangles indicate the 
patches of scintillation counters of the MIA underground array (64 counters/patch). Five crossed 
circles indicate tracking Cherenkov telescopes. The D.I.C.E. detectors are not used in this analysis. 
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Figure 3: Shadow of the Moon as detected by CASA-MIA. The contour plot shows the deficit 
in the number of detected cosmic ray events as a function of the angle from the Moon center. 
The axes are defined by the equatorial coordinates of the Moon, with right ascension along the 
horizontal axis and declination along the vertical axis. The contours correspond to successive one 
standard deviation steps. A Gaussian smoothing factor of 0.5° has been applied to the data. 
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Figure 4: Angular resolution of CASA-MIA as a function of the number of alerted CASA stations. 
The resolution has been estimated by three different techniques: split-array method (solid points), 
Cherenkov telescope array coincident events (open triangles), and Moon shadow (open squares). 
The dashed curve indicates a simple parametric fit to the data. 
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Figure 5: Histograms of the number of muons per event detected by CASA-MIA for cosmic ray 
data events (solid) and simulated gamma-ray events (dashed). The gamma-ray simulation uses 
an input energy spectrum with the same power law spectral index as the cosmic ray data. 
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Figure 6: Distributions of the relative muon content, r M , for cosmic ray data events (solid) 
and simulated gamma-ray events (dashed). The quantity r M is defined in the text. Events with 
zero muons are assigned a value of r^ = —3.5. The gamma-ray simulation uses an input energy 
spectrum with the same power law spectral index as the cosmic ray data. 
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Figure 7: Scans in right ascension for a band of constant declination centered on Cygnus X-3, 
at a median energy of 115 TeV. The data points correspond to the numbers of events observed 
in each bin; the dashed histogram is the expected background level. Top plot shows the all-data 
sample; bottom plot shows the muon-poor sample. Note that the scale on the vertical axis has 
been highly expanded and zero-suppressed. 
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Figure 8: Scans in right ascension for a band of constant declination centered on Hercules X-l, 
at a median energy of 115 TeV. The data points correspond to the numbers of events observed 
in each bin; the dashed histogram is the expected background level. Top plot shows the all-data 
sample; bottom plot shows the muon-poor sample. Note that the scale on the vertical axis has 
been highly expanded and zero-suppressed. 
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Figure 9: Distribution of transit significances for Cygnus X-3 for the all-data (top) and muon- 
poor (bottom) samples. The curves are unit-width Gaussian distributions with zero mean. For 
comparison, the fitted mean (standard deviation) of the all-data sample is —0.040 ±0.031 (1.011 ± 
0.021) and the fitted mean (standard deviation) of the muon-poor sample is —0.011 ± 0.032 
(0.992 ±0.022). 
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Figure 10: Distribution of transit significances for Hercules X-3 for the all-data (top) and muon- 
poor (bottom) samples. The curves are unit-width Gaussian distributions with zero mean. For 
comparison, the fitted mean (standard deviation) of the all-data sample is —0.032 ±0.031 (0.973 ± 
0.020) and the fitted mean (standard deviation) of the muon-poor sample is 0.030 ±0.033 (0.976 ± 
0.022). 
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Figure 11: Phase distributions for events from Cygnus X-3 based on the 4.8 hr X-ray periodicity. 
Data points show the numbers of on-source events in each phase bin for the all-data (top) and 
muon-poor (bottom) samples. Histograms show the expected background levels. 
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Figure 12: Flux limits reported between 1990 and 1995 on the steady emission of particles from 
Cygnus X-3. The squares (5) represent the results of this work. Open squares indicate limits on 
the emission of any neutral particle that creates air showers. Filled squares indicate limits on the 



emission of gamma-rays. The circles (1-4) represent results from other experiments: 1. Tibet [92], 
2. HEGRA H, 3. CYGNUS @, and 4. EAS-TOP @. The dashed curve is the approximate 
power law fit to early results (reproduced from Figure |l]). 
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Figure 13: Limits reported between 1990 and 1995 on the fractional excess of gamma-rays from 
Cygnus X-3 relative to the cosmic ray background. The squares represent the results of this work. 
Open squares indicate limits on the emission of any neutral particle that creates air showers. 
Filled squares indicate limits on the emission of gamma-rays. The circles represent results from 
other experiments: Tibet @, HEGRA ||, CYGNUS |j§, and EAS-TOP @. The dashed 
lines indicate the range of fractional excess values corresponding to the fluxes reported by earlier 
experiments. 
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Figure 14: Flux limits reported between 1990 and 1995 on the steady emission of particles from 
Hercules X-l. The squares (5) represent the results of this work. Open squares indicate limits 
on the emission of any neutral particle that creates air showers. Filled squares indicate limits on 
the emission of gamma-rays. The circles (1-4) represent results from other experiments: 1. Tibet 
||, 2. HEGRA m, 3 - CYGNUS H, and 4. EAS-TOP @. 
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Figure 15: Daily event totals in 1994 for showers recorded by CASA-MIA from the direction of 
Hercules X-l for the all-data (top) and muon-poor (bottom) samples. The dashed lines indicate 
the expected numbers of events for gamma-ray fluxes comparable to those seen in 1986 |j(J . The 
gaps in the distribution indicate those periods in which the experiment recorded no usable data 
from Hercules X-l. 
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